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How do fermions behave on a random lattice ? 

C. J. Griffin and T. D. Kieu* a 

a School of Physics, University of Melbourne 
Parkville, Victoria, Australia, 3052 

Comparing random lattice, naive and Wilson fermions in two dimensional abelian background gauge field, we 
show that the doublers suppressed in the free field case are revived for random lattices in the continuum limit 
unless gauge interactions are implemented in a non-invariant way. 



1. Introduction 

The doubling problem of lattice fermions is in- 
evitable according to the Nielsen Ninomiya no- 
go theorem[[l] if the free-field action satisfies the 
conditions of reflection positivity, locality, global 
axial symmetry, and translational invariance at 
a fixed scale. An obvious resolution of the dou- 
bling problem is thus to relax one of those con- 
ditions to obtain, in the order listed above, non- 
hermitian[|, non-local||], Wilson[Q, or random- 
lattice p-p fermion formulations. These formu- 
lations are all free of doublers when there are 
no interactions or when the interactions are of a 
non-gauge nature |J: the extra poles in the prop- 
agators are removed as the lattice spacing a de- 
creases, leaving a single fermion mode in the con- 
tinuum limit. 

Gauge interactions behave very differently on 
account of a unique and special property. Lo- 
cal gauge invariance imposes severe constraints 
on the theory, expressed mathematically in the 
Ward-Takahashi identities. In particular, the 
fermion-gauge vertex is related to the free inverse 
propagator, 



V» ~ d M G (p), 



(1) 



giving the interaction vertices mode dependency. 
This different coupling strength of doublers to 
gauge fields has been shown to revive these 
modes in loop diagrams, even though they are 
suppressed at the free-field level, in studies of 
some non-local and non-hermitian formula- 
tions [jll 12 1 . For this reason, we investigate the 
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issue of fermion doubling on random lattices with 
gauge interactions |l3j. 

In the random lattice approach, suitable quan- 
tities are measured on a random lattice then av- 
eraged, cither qucnchcdly or anncalcdly, over an 
ensemble of lattices. Apart from the extra work 
involved in generating an ensemble of random 
lattices, this approach better approximates the 
scale-free rotational and translational symmetry 
of the continuum than that of regular lattices. 
Thus, the scaling region near the continuum limit 
may be more easily reached on random lattices 
than on regular lattices of the same size. More 
relevant to this discussion, since there is no fixed 
Brillouin zone, there need be no extra poles of the 
propagator. Even if extra poles do exist, the one- 
to-one correspondence between propagator poles 
in momentum space and zero modes is not neces- 
sarily valid since plane waves are no longer eigen- 
states of the Dirac operator ||. 

This expectation of no doubling on random lat- 
tices has been confirmed in various studies of free- 
field theory in both two and four dimensions 
It has similarly been shown that random lattice 
theories with four-point interactions are also dou- 
bler free|]]. 

2. Our Approach 

We wish to compute the fermion determinant 
of abelian background gauge theory on a two di- 
mensional Euclidean random lattice which we ap- 
proximate with 

- InDet (GaGo 1 ) = Tr [G^G - l] - 
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Figure 1. Typical lattice structure 



-Tr [(G^Go-l^+OG? 4 ), (2) 

where G^ 1 is the fermion propagator in the back- 
ground gauge field A^. Background fields are the 
most appropriate for our purposes, providing a 
clean signal which increases with the number of 
fermion flavours contributing to the internal lines. 
Comparing with identical calculations for naive 
and Wilson fermions on square lattices, which 
are known to be four-fold doubling and doubler- 
free respectively, clarifies the continuum limit be- 
haviour of our random lattices. 

Our lattice is constructed from a triangulated 
array of N square lattice vertices by a random se- 
quence of Alexander "flip" moves, supplemented 
by further constraints which force the lattice to 
stay locally flat throughout the flipping proce- 
dure. We chose to randomise the lattice with 6N 
successful flips. See figure [j]. This fixed vertex 
construction has several advantages over the lat- 



tices of Construction is O(N), compared with 
0(N 3 ), and the resulting lattice has a fixed vertex 
spacing (a fixed size), so measured quantities do 
not need to be adjusted for different link-lengths. 
The trade off is that we have introduced a lattice 
dependent internal scale, s, which will need to be 
dealt with. 

The fermion derivative is chosen such that it 
reduces to the naive result on lattices of regular 
arrangements of links. It is constructed by av- 
eraging the contributions of pairs of consecutive 
links (fc, I) around each vertex to the derivative 
in lattice framing co-ordinates, replacing the con- 
tinuum T.^l^^i. 30 ) w i tn 

- V tt—, x [I ^x+k - k y x +i + (k-l) * x ] 

i— l...n 

Gauge interactions are introduced in the usual 
gauge-invariant manner using the link variables 
Ux.x+i = exp(i<7 J ; A(x) . Ax). An alternative defi- 
nition U x , x +i = exp(igl . A(x + l/2)) which allows 
gauge invariance to be explicitly broken on the 
lattice, is also considered. The resulting action is 
hermitian, local and axially-symmetric. 

Measurements are made in a background gauge 
field 

with 

fixed physical quantities raa -1 = 0.1 (length) -1 , 
a 2 N = 64 (length) 2 , and a~ 2 £ = 0.05 (length)" 2 , 
for a = {1.0,0.5,0.3333,0.25} 

3. Results 

We first compute a quantity derived from the 
free propagator 

/(0 = Tr 7 1 I (l + 7o )G -W) x 

S^xo-x'o+O, (4) 

evaluating the average zero momentum real par- 
ticle propagator projected along the xq direction. 
Figure ^ summarises the calculation, clearly iden- 
tifying the doubler suppression of free fermions 
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Figure 2. Fermion propagation, /(£). 



Figure 3. Variation of — InDet (GaGq w i tn a 



on a random lattice in agreement with jf). In- 
deed, apart from some minor small distance fluc- 
tuations of the order of the internal scale, the ran- 
dom lattice result matches the continuum com- 
pletely; the normalisation is reproduced exactly, 
and masses do not need to be tuned. 

The calculation of the approximation in equa- 
tion ^| is complicated a little by its sensitivity 
to the internal scale of the lattice. Ideally we 
would like to ignore this effect as a small correc- 
tion, or have the internal scale match the lattice 
spacing (s) ~ a. In practice, (s) ~ 1.3a, the 
contribution is significant, and (s) gets larger as 
the degree of randomisation is increased. Using 
an ensemble of lattices, it is possible to extrapo- 
late to (s) = a, leading to a mean value of this 
extrapolation, and uncertainties associated with 
the spread of geometrically different lattices that 
have the same (s). See figure ||| for the results. 
The naive case approaches the continuum limit 
quadratically with a and the Wilson approaches 
~ 4 of the same result linearly, as expected. The 



same graph also indicates the random lattice re- 
sults; the number of lattice configurations used 
in the extrapolation is displayed next to each 
point. The lattice gauge invariant calculation 
is clearly more like the naive fermion than the 
Wilson. Moreover, had extrapolation in (s) not 
been performed, the result would have been even 
larger. With gauge invariance broken, the con- 
verse is clearly seen; the result is certainly more 
like Wilson than naive, and in this case the ex- 
trapolation procedure has forced an increase. An- 
other observation is that the gauge non-invariant 
calculation approaches the continuum result more 
rapidly than either Wilson or naive formulations. 

4. Discussion 

It is clear from our results that there are dou- 
blers on random lattices when gauge invariance 
is maintained at finite lattice spacing, since the 
extrapolated determinant is comparable to that 
of naive fermions. It can also be seen that the 



4 



doubling can be avoided if one gives up gauge 
invariance (but needs and hopes to recover it in 
the continuum limit, as is in the case we studied 
above). 

In all cases, the lattice fermion actions are in- 
variant under the global axial transformations. 
When there are doublers on random lattices, the 
axial anomalies are canceled in the usual man- 
ner among opposite-chirality species. When there 
is no doubling in the gauge non-invariant for- 
mulation, the conserved lattice current being the 
Noether current of axial symmetry is, of course, 
not gauge invariant. Thus it cannot be identi- 
fied with the continuum axial current which is 
invariant. It should be, instead, identified with 
a combination of the continuum current and a 
gauge-noninvariant term, whose divergence gives 
us the axial anomalies, 

= JZL^to+a^M*)- (5) 

We believe that the results obtained here are 
also applicable to other kinds of random lattices. 
Our doubling conclusion for random lattices is not 
plainly disappointing but also points to some se- 
rious implications. 

We have extended the lattice no-go theorem 
and at the same time emphasised the importance 
of gauge invariance in the phenomenon of lattice 
fermion doubling. 

The failure of random lattices to accommodate 
chiral fermions could either undermine the point 
of view that at the Planck scale or higher the 
structure of spacetime is that of randomness; or, 
taken with other complete failures in dealing with 
chiral fermions, could be a hint that our under- 
standing of chiral gauge theories is incomplete. 
And, correspondingly, the quantisation of those 
theories is in need of further studies. One of us 
has been pursuing this latter path]l4]|. 
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